home *** CD-ROM | disk | FTP | other *** search
/ Cream of the Crop 1 / Cream of the Crop 1.iso / PROGRAM / NRPAS13.ARJ / TPTEST.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  586b  |  22 lines
  1. PROCEDURE tptest(data1,data2: glnparray; n: integer; VAR t,prob: real);
  2. (* Programs using routine TPTEST must define type
  3. TYPE
  4.    glnparray := ARRAY [1..np];
  5. in the main routine, with np >= n. *)
  6. VAR
  7.    j: integer;
  8.    var2,var1,sd,df,cov,ave2,ave1: real;
  9. BEGIN
  10.    avevar(data1,n,ave1,var1);
  11.    avevar(data2,n,ave2,var2);
  12.    cov := 0.0;
  13.    FOR j := 1 TO n DO BEGIN
  14.       cov := cov+(data1[j]-ave1)*(data2[j]-ave2);
  15.    END;
  16.    df := n-1;
  17.    cov := cov/df;
  18.    sd := sqrt((var1+var2-2.0*cov)/n);
  19.    t := (ave1-ave2)/sd;
  20.    prob := betai(0.5*df,0.5,df/(df+sqr(t)))
  21. END;
  22.